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Abstract 

We study Matrix Element corrections as implemented in four popular 
event generators for hadron collisions. We compare PYTHIA, HERWIG, AlpGen 
^ ! and SHERPA in the Z/7* inclusive production at LHC. PYTHIA and HERWIG are 

00 \ able to correct the first emission from the shower taking the Matrix Element 

calculation for one additional parton into account. SHERPA and AlpGen are 
. able to take into account Matrix Element corrections not only for one, but 

\ rather for several hard emissions from the incoming partons. This can be 

^) ' done at the price of introducing a separation cut to distinguish a Matrix 

Q!^ ■ Element and a Parton Shower populated regions. In this paper we check 

the effect of Matrix Element corrections in PYTHIA and HERWIG and we check 
that results from these two generators are consistent. Then we turn to SHERPA 
^ ' and AlpGen, that implement two different methods to match Matrix Element 

^ . calculations and Parton Shower. If we constraint them so that no more than 

one parton can emerge from the Matrix Element calculations they should 
both give results similar to PYTHIA and HERWIG. In other words PYTHIA and 
HERWIG provide us with the correct reference to spot possible issues with 
the matching prescriptions implemented in SHERPA and AlpGen. We also 
check to what extent the dependency on the Matrix Element - Parton Shower 
separation cut is canceled in these two generators. 
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1 Introduction 



The description of the QCD radiation pattern that accompanies the partons in- 
volved in a hard scattering process can be done at different precision levels. Most 
event generators attach the QCD radiation to the partons using the Parton Shower 
technique. While this approach provides a good description of many low pr ob- 
servables it may fail in efficiently filling the phase space for hard radiation. This 
limitation is connected to the Parton Shower being a collinear approximation of the 
decription of parton splittings. 

A way to improve the description of the QCD radiation pattern is to augment 
the Parton Shower with information coming from the exact matrix element calcula- 
tion. The generators studied in this work follow this approach, but with significant 
differences. PYTHIA [T] and HERWIG [21 E] are able to modify the shower in such 
a way that the hardest emission is described using the exact Matrix Element cal- 
culation for one additional real emission. AlpGen [1] and SHERPA [5] are able to 
take Matrix Element corrections for several (not only one) hard emissions. The 
main idea in both generators is that configurations in which the emitted partons' 
Pt is below a certain threshold are described with a pure Parton Shower approach, 
while configurations in which n partons are above the threshold are described with 
the n-real emissions Matrix Element. So, at the price of introducing an arbitrary 
threshold that is not present in PYTHIA or HERWIG, AlpGen and SHERPA should be 
able to describe the emission of several hard partons with the corresponding Ma- 
trix Element calculation. The dependency on the arbitrary threshold has to be as 
limited as possible. 

In this paper we compare these different approaches for the case of the inclusive 
Z/7* production at LHC. We first check the effect of Matrix Element corrections in 
PYTHIA and HERWIG. We discuss the differences in the implementation and compare 
results from these two generators. Then, we discuss the differences in the merging 
prescriptions implemented in SHERPA and AlpGen and we present a sort of a consis- 
tency test about these prescriptions. The main idea behind this test is that if we 
allow the Matrix Element generators in SHERPA and AlpGen to emit not more than 
one additional parton we should recover results from HERWIG and PYTHIA, because 
we are using the same Matrix Element content. We also check to what extent the 
dependency on the threshold value in AlpGen and SHERPA is canceled. 

2 Improving the description of the QCD 
radiation pattern using Matrix Element 
Calculations 

In this section we briefly describe the implementation of the Matrix Element cor- 
rections in the generators studied in this work. We group the PYTHIA and HERWIG 
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implementations under the name of "Parton Shower reweighting" , while SHERPA 
and AlpGen belong to the cathegory of "Matrix Element - Parton Shower merging" . 

2.1 Parton Shower reweighting 

In this section we describe the implementation of Matrix Element corrections in 
the Initial State Parton Shower as implemented in PYTHIA. We will enlight the 
differences between HERWIG and PYTHIA at the end of this section. 

Matrix element corrections for inclusive Zj^* production in PYTHIA modifies 
the first emission from the PS so that the effect of the first order correction is 
reproduced. The lowest order graph contributing to Zj^* production is shown in 
Fig. [TJ The processes that contribute to the first order correction are qq — > Zj^*g^ 
represented by the graphs in Fig. [2l and qg Zj^*q^ represented by the graphs in 
Fig. El 

The correction comes as a re-weighting factor of the parton shower. The initial 
state parton shower is weighted with two factors, one to reproduce the matrix 
elements of Fig. [2], Wqq^zj'-^'g-, and one to reproduce the matrix element of Fig. [HI 
Wqg^zh'q [6]. The PS emission closest to the hard qq process is the one that gets 
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Figure 1: Lowest order contribution to the Zj^* production. 




Figure 2: Graphs contributing to the qq Zj^Q process. 

the correction. In order to correct the PS we need to classify this emission as either 
[qq — ^ Z/'-f*g) -like or {qg Z/'~f*q) -like. A PS branching like the one depicted in 
Fig. m is considered {qq —>■ Z/'~f*g) -like, while a branching like the one in Fig. is 
considered {qg Z/7*g)-like. 
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(a) (b) 
Figure 3: Graphs contributing to the qg Z/j*q. 




Figure 4: The circled PS branching is considered of type qq — > Z/'j*g. 




Figure 5: The circled PS branching is considered of type qg — > Z/'~f*q. 



The ratio of the ME and PS differential cross sections, classified as described 



above, leads to the following expressions for 



qq^z/rg and Wgg^z/rq 0: 



w, 



_ P + + 2m|s 
qg^z/rq = (s_rn|)2 + m|' 



(1) 
(2) 



where s, t and u are the Mandelstam variables. It can be shown that 



< W, 



qq^Z/-y*g < 1, 



and 



1 < W, 



.„ < 3. 



(3) 
(4) 
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This means that the Parton Shower alone would overestimate the qq Z/'-)*g 
contribution and underestimate the qg —>■ Z/'j*q contribution. 

As described in |6] the correction is performed in two steps: first the starting 
scale for the shower is raised to the hadronic center of mass energy, so that any hard 
emission from the shower is kinematically possible; then the first splitting from the 
shower is classified as {qq — > Z/'j*g)-\ike or {qg Z/'y*q)-\ike and the splitting 
probability is modified according to the corresponding W function. 

Matrix element corrections are implemented in HERWIG in a similar way, but 
with two important differences: 

• in HERWIG the PS is angular ordered, which means that the early emissions 
are soft large angle gluons; thus, the emission that deserves the corrections is 
not the first like in PYTHIA, as pointed in jZj; 

• in HERWIG the PS cannot fill the phase space for values above the hard scale 
of the process; in this respect it is similar to PYTHIA with starting scale set to 
Mz. 

For these reasons a two-step approach is used in HERWIG to implement ME correc- 
tions [8]. In the phase space region covered by the PS, corrections are applied as in 
PYTHIA, the only difference being that not the first emission, but rather the hardest 
emission so far during the shower evolution gets the correction [7] . These correc- 
tions are referred to as "soft ME corrections". In the "dead zones", that are left 
completely uncovered by the PS, the exact ME for one additional parton is used, 
with subsequent PS. These corrections are referred to as "hard ME corrections". 

2.2 Matrix Element - Parton Shower merging 

SHERPA and AlpGen implement matrix element corrections in a different ways. They 
both subdivide the phase space in two regions, one for jet production which is filled 
by the Matrix Element calculation, and one for the jet evolution, which is filled by 
the Parton Shower. The criterion that is used to distinguish the two regions is a jet 
measure based cut. 

The procedure implemented in SHERPA is known as CKKW merging. It was first 
proposed for lepton collisions in [9j , then it was extended to hadron collisions in [TUj . 
In a nutshell, the prescription foresees a preliminary step in which cross sections 
are calculated for up to n additional partons in the final state; ME level events 
are produced according to the calculated cross sections and they are weighted with 
a Sudakov weight that makes these states exclusive. In other words the Sudakov 
weight accounts for the probability that, given an n parton final state, no further 
hard emission is done by the Parton Shower; the reason for this is that any additional 
hard emission is going to be simulated by the n + 1 Matrix Element. 

More extensively, the whole procedure foresees the following steps: 
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1. ME cross sections (j„.j are calculated for each parton multiplicity n and for 
each different combination i of partons that contributes to multiplicity n. A 
cutoff Hcut on the separation of partons is applied to avoid divergences. A 
fixed as^^^ is used. 

2. One among configurations n, i is selected with probability Pn,i = cr„_j/ Ylim j ^"i, 

3. Parton momenta are generated according to the corresponding matrix element 
squared. 

4. The scales at which the splittings happened are reconstructed: this is achieved 
through a k± clustering of the partons emerging from the ME. The clustering 
is stopped when the core 2 — >• 2 process is found. This leads to a series of 
n — 2 clusterings with associated values of the k± distance j/2---2/n- Once the 
values yi = tf/s are known we can finally calculate the ME event weight, that 
comes in two factors: 

- an as correction: for each clustering i an as correction as{ti)/a^^^ is 
applied; 

- a Sudakov form factor correction is applied, to make states exclusive. 

5. Events are accepted or rejected according to their weight. 

6. The accepted events are showered with a veto on the emission above ycut- 

The matching prescription implemented in AlpGen is known as MLM [TT]. It 
is similar in the motivations to the CKKW, but different in the implementation. 
Basically, instead of calculating Sudakov weights analytically, the event weighting 
is performed numerically. As in CKKW, ME cross sections are calculated up to the 
maximum parton multiplicity that the user wants in the final state; a minimum px 
cut for final state partons is used to cutoff ME divergences and a fixed as is used. 
As in CKKW a "PS history" is reconstructed and a splitting sequence is identified, 
with corresponding scales; an as correction is applied as in the CKKW. From this 
point the two prescriptions become different. In the MLM approach a conventional 
PS program is used (PYTHIA or HERWIG) and ME partons are showered without any 
constraint. The parton collection that results from the PS step is clustered using 
a jet algorithm (a cone in the AlpGen implementation, but also other options have 
been investigated, e.g. a k± algorithm is used in MADGRAPH/MADEVENT) ; the resulting 
jets are matched in angle to the ME partons and only those events in which all the 
jets match the ME partons without any extra unmatched jets are retained (for the 
maximum ME parton multiplicity additional jets, softer than the matched ones, are 
allowed) . 

This procedure tries to reproduce in one go the effect that in the CKKW is 
achieved in two steps: the Sudakov re-weighting and the vetoed shower. Indeed 
the rejection of events with additional jets should, at the same time, reject ME 
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configuration in a similar way as the Sudakov weight does and prevent additional 
emission from the shower, thus reproducing the effect of the CKKW PS veto. 

The MLM prescription is really convenient because it does not require modifi- 
cations in the PS program. It just requires a veto routine to kill events not fulfilling 
the matching criteria. 

While the CKKW prescription contains one parameter (the qcut of the internal 
clustering algorithm), in the MLM the user has to choose different parameters. The 
cone algorithm used for the matching has three parameters, namely the minimum 
jet pt the cone radius i?, and the jet maximum pseudorapidity Tj. The minimum 
Pt used in the cone clustering (pTmL) '^'^^ ^^e same as the minimum pt used 
in the ME step to cutoff divergences ipTmin)'- usually it is recommended to have 
PTmin >PTmiri'i ^^is is needed because events that are below the cut at the ME level 
could fall above after the PS. For this reason a process dependent tuning for PTmin 
with respect to PTmIn needed; for the Z/Y + jets production AlpGen authors 
recommend to choose the jet finder minimum pt to be 5 GeV higher than the ME 
minimum pt- 



3 Generator configuration and analysis 



Plots shown in this paper were obtained with an analysis coded in the Rivet fTIi [TH] 
MC Validation/Tuning framework. Rivet comes with an interface library, called 
AGILe, which provides interface to many event generators. The user can run the 
event generator through the AGILe interface and analyse the HepMC event record 
exploiting Rivet analysis classes. While Rivet can be used as an analysis framework 
for general studies on Monte Carlo event generators, as we did in this paper, its 
main feature is that it comes as a very straightforward tool to make comparisons 
to data. 

In order to gain a more detailed understanding of the hard event simulation 
all the generators used in this work were run switching off the underlying event 
simulation. The analysis has been carried out at parton level. Also, the emission 
of final state photons from the leptons from the Z boson decay has been switched 
off. In the analysis, the lepton pair from the Z boson decay was required to have 
an invariant mass between 66 and 116 GeV. Jets were reconstructed with the longi- 
tudinally invariant kj_ algorithm [H], as implemented in the Fast Jet [12] package. 
The pseudo-radius parameter in the k±_ algorithm was set to 0.4 and the minimum 
Pt for jets was set to 30 GeV. 
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4 Matrix element corrections in PYTHIA and HERWIG 



4.1 Z boson transverse momentum 

The distribution for the lepton pair for inclusive Zj^*' production in PYTHIA is 
shown in Fig. [HI Only the electron decay channel has been selected. 




Figure 6: Pt spectrum for the lepton pair in PYTHIA for three different settings of 
the ISR: with ME corrections, without ME corrections and with the starting scale 
of the shower set to y^, without ME corrections and with the starting scale of 
the shower set to Mz- The relative difference with respect to the curve with ME 
corrections is shown in the lower plot. 

The three curves correspond to three different configurations: one is with ME 
corrections activated, the other two are obtained without matrix element correc- 
tions, but with different starting scales for the shower: the total hadronic center of 
mass energy and the invariant mass of the lepton pair respectively. 

PYTHIA implements a virtuality ordered parton shower. The starting scale of 
the shower marks the maximum allowed virtuality in the shower evolution. If the 
starting scale is set to Mz the hardest parton transverse momentum cannot ex- 
ceed Mzi thus also the Z cannot exceed 90 GeV approximately, as shown in 
Fig. [HI When the starting scale is raised, the spectrum gets harder. When ME 
corrections are activated the spectrum gets even harder. The reason why the ME 
corrected spectrum is harder than the uncorrected one can be explained consider- 
ing the relative amount of the two corrections at the LHC The graph gq Z/'y*q 
contributes more than qq Z/'~f*g, because q is not a valence quark at LHC. We 
recall that Eq. (jH) states that Wgy^z/^/'q > 1, meaning that the first emission from 
the ME-corrected shower is done with a splitting probability higher than that of 
the uncorrected shower. Since the PS emission is always ordered, a higher splitting 
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probability means that the probabihty for emitting a harder parton is higher than 
that of uncorrected shower. This explains why the corrected spectrum is harder 
than the uncorrected one. We also notice in passing that at the Tevatron pp col- 
hder [16l ITTl HHj exactly the opposite holds: in that case qq Z/ ^9 dominates 
because g is a valence quark; since the correction is < 1 for this graph (Eq. ([3])) the 
ME corrected result is softer than the uncorrected one, as shown in [6]. 

One might expect that ME corrections should change the shape of spectra only 
at high pt'i low pt region should be well described by the parton shower alone. 
Actually the Z px spectrum is altered by ME corrections all over the px range, as 
shown in Fig. [HI The three distributions are normalized, but the difference at low 
Pt is not only due to normalization; a change in shape is also present; such a change 
is testified by the relative difference plot, that does not flatten as px approaches 
zero. The reason for the change in shape at low px is that ME corrections change 
the Sudakov form factor used in the shower [TO]. 

The Z Pt spectrum as obtained in HERWIG is shown in Fig. [71 Both the ME- 
corrected and the uncorrected spectra are shown. While ME corrections in PYTHIA 
change the whole shape of the distribution, also at low pt, the change in shape at 
low Pt in HERWIG is small. 

A comparison between ME corrected Z pt distribution in PYTHIA and HERWIG 
is shown in Fig. [HI The agreement is very good all over the pt spectrum. 
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Figure 8: ME corrected Z pt distribution in PYTHIA and HERWIG. 



4.2 Differential jet rates 

To test how the phase space available for QCD radiation is filled we looked at 
differential jet rates. The differential jet rate n n — 1 is the distribution of the 
n —>■ n — 1 transition value, Qn^n-i, which is the value of the resolution parameter 
dcut (in an exclusive k±_ algorithm [T3]) for which an n jet event turns into an ra — 1 jet 
event. To understand what a differential jet rate is let's consider a simple example 
with three particles in the final state. Let the particles be 1, 2, 3. Let di^B be the 
k± distance of particle i from the beam line and dij be the k± distance between 
particles i and j. Suppose that the sequence of ordered distances looks like this: 



dl,B < d2,B < rfl,2 < rfs.B < C?l,3 < C?2,: 



(5) 



In an exclusive calculation the first jet that would be recombined with the beam 
line is particle 1, thus the 3^2 transition value is (^3^2 = di^B- Then particle 2 is 
the next one to be recombined with the beam, thus (^2^1 = '^2,_b- At this point, if 
the dcut is raised to be at least (i2,s only d^^B survives in Eq. ([5]), thus Qi-^o = d^ b- 
Let's now consider another example sequence: 



dl,B < C?l,3 < C?l,2 < C?2,3 < d2,B < ds^B- 



(6) 



As before, Q3^2 = c^i.b- If we raise dcut to at least di^B particle 1 gets clustered 
with the beam and the new sequence looks like 



(^2.3 < d2,B < d. 



(7) 



so particles 2 and 3 are going to be clustered in the next step. Thus if dcut is set 
to be at least equal to ^2,3 particles 2 and 3 are clustered, thus passing from 2 to 
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1 jet. Thus Q2^i = c?2,3 After they are clustered the sequence will have one only 
element, namely d(2fi),B] thus Qi^o = '^(2,3),b- 

Let's now see how the differential jet rates in PYTHIA look for the three setting 
of ME corrections. The 1 — > differential jet rate is shown in Fig. O The main 




differences among the three settings are in the region for high values of Qi^q. We 
see that the sample without ME corrections and with low starting shower scale is 
the one that dies first. This means that it is unable to fill the phase space for 
hard parton emission, which is responsible for filling the rightmost part of the plot. 
Regarding the other two settings we observe that they are similar, but while the 
sample without ME corrections tends to fill the region below 1, the one with ME 
corrections fills the region above 1 more, thus allowing for more radiation to be 
emitted. 



5 SHERPA and AlpGen compared to PYTHIA 

As mentioned in the introduction, we made a consistency check of the matching 
prescriptions implemented in SHERPA and AlpGen. The test consists in the compar- 
isons of the observables obtained with ME corrected PYTHIA with those obtained 
with AlpGen and SHERPA when at most one additional parton is allowed to come 
from the ME calculation. This approach has been already explored in to test 
various matching prescriptions in e^e~ collisions. 
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5.1 SHERPA 



The total cross section for inclusive Zj^* is calculated in PYTHIA at Leading Order 
accuracy, i.e. it is calculated from the process qq — > e+e~; the subsequent PS, either 
corrected or uncorrected, happens with probability 1, so it cannot modify the cross 
section. In SHERPA the cross section comes as a sum over the different selected final 
state parton multiplicities. 



CXSHERPA = ^^Oi <Wi> (8) 
j=0 



where Oi is the cross section for % additional partons in the final state and < Wi> 
is the average Sudakov weight for that configuration. 

Table [U shows the cross section values as obtained in PYTHIA and in SHERPA for 
four different values of the resolution cut q^ut that steers the separation between the 
ME and the PS regions. 



SHERPA 




cfi [pb] 


<Wi> 


Total o [pb] 


g,„i=10GeV 


jet 


838.9 


0.7489 


1383 




1 jet 


998.7 


0.7559 




gent=20GeV 


jet 


1059.5 


0.9301 


1405 




1 jet 


484.6 


0.8657 




gent=40GeV 


jet 


1271.2 


0.9992 


1434 




1 jet 


177.2 


0.9267 




ge«t=500GeV 


jet 


1926.6 


0.7540 


1453 




1 jet 


0.038 


0.9802 





PYTHIA Total o [pb] 



inclusive 1528 



Table 1: Cross sections for SHERPA and PYTHIA. 



The difference in the total cross section with respect to PYTHIA is up to about 
10%, for the sample with the lowest value of qcut- Cross sections for both SHERPA and 
PYTHIA are formally LO; some differences are due to the qcut dependency mainly. 
If one takes a very high value for q^uti this makes SHERPA get closer and closer 
to PYTHIA. In fact, as q^ut is increased the contribution to the total cross section 
from the configuration with one additional parton vanishes, thus leaving the leading 
order contribution alone, which is the only one considered in PYTHIA for the cross 
section calculation. In summary the dependency on qcut in SHERPA is of the order 
of 5-10%. Even when higher order emissions are completely removed there is a 5% 
discrepancy with PYTHIA. This presumably happens due to differences in the choice 
of scales in the two generators. In particular in SHERPA the renormalization scale is 
set to the qcut value. 
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Figure 10: Pt spectrum (a) and rj distribution (b) for the lepton pair in PYTHIA and 
SHERPA. The latter has been run with at most one additional parton from the ME; 
three different values for the separation cut between ME and PS regions have been 
used: qcut= 10, 20, 40 GeV. 

The Pt spectrum for the lepton pair in SHERPA and PYTHIA is shown in Fig. [10] 
(a). As mentioned above, SHERPA has been run such that only one additional parton 
can be emitted from the matrix element. 

SHERPA has been run with three different values of the parameter qcut that gov- 
erns the separation between the phase space region filled by the ME and the region 
filled by the PS. The values used were qcut= 10, 20, 40 GeV. As described in Sec. 12.21 
the region above Qcut is filled by a modified ME, in which Sudakov form factors are 
attached to the ME, while the region below the cut is filled by a vetoed PS. SHERPA 
appears to closely follow the PYTHIA spectrum. Some discrepancies are observed in 



13 




Figure 11: pr distributions for the positive lepton. 



the first few bins. These are most probably due to a different tuning for the pri- 
mordial ptH distribution of partons in the protons. We notice in particular that the 
high pt tail, that is sensitive to additional hard emission, appears to be correctly 
reproduced. We also observe a small dependence on the value of qcut- the three 
curves for SHERPA agree within few percent. 

The pseudorapidity rj distribution for the lepton pair is shown in Fig. [TU] (b). 
The distributions from SHERPA well agree with PYTHIA in the central region, where 
the difference is within 10-15%; in the tails of the distribution difference is more 
evident. 

The Pt distribution for the positive lepton is shown in Fig. [TTl In this case too 
SHERPA agrees with PYTHIA, without strong dependency on the resolution cut qcut- 

Concerning the QCD observables we looked at differential jet rates. The distri- 
bution for the rate 1 ^ in PYTHIA and in SHERPA is shown in Fig. [T21 For SHERPA 
we tried three different values for the matching parameter qcut- 

10, 20, 40 GeV. The 

vertical lines shown in the plots indicate the position of the resolution parameter 
qcut- In SHERPA, the phase space above qcut is filled by the Matrix Element, while 
the region below qcut is filled by the vetoed shower. 

SHERPA agrees well with PYTHIA. The transition between the ME- and PS- 
populated regions is quite smooth for all the three values used for qcut- The relative 
difference with respect to PYTHIA is at most 20%, in the vicinities of the qcut region. 
The tail for very low values of Qi^o shows some differences, that can be due to 
slightly different settings for the primordial k±. 



^The primordial pT distribution of partons in protons, often referred to as primordial fcj^ is the 
transverse momentum distribution of partons in the hadrons entering the collision. The primordial 
k± models a soft non-perturbative effect. 
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Figure 12: Differential jet rate 1 ^ in PYTHIA and SHERPA. For SHERPA we used 
three values for the matching parameter qcut'- 10, 20, 40 GeV. Relative differences 
are calculated with respect to PYTHIA. 

5.2 AlpGen plus PYTHIA 

We made the same test with AlpGen, looking for differences with respect to ME 
corrected PYTHIA. We considered the contribution from up to one additional parton 
from the matrix element, and we used PYTHIA to shower the ME events generated 
by AlpGen. 

Table [2] summarizes the total cross section for AlpGen with up to one additional 
parton from the ME. Results for three values of the ME cutoff for the generation of 
the additional jet are shown. The difference with respect to PYTHIA is 5% at most. 

The Pt spectrum and the t] distribution for PYTHIA (with and without ME 
corrections) and AlpGen are shown in Fig. [121 Concerning PYTHIA, the sample 



AlpGen 




(Ti [pb] 


Total a [pb] 


ME cutoff=10GeV 


jet 


1092 


1609 




1 jet 


517 




ME cutoff=20GeV 


jet 


1303 


1594 




1 jet 


291 




ME cutoff=40GeV 


jet 


1452 


1558 




1 jet 


106 




PYTHIA 






Total a [pb] 


inclusive 






1528 



Table 2: Cross sections for AlpGen and PYTHIA. 
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Figure 13: px spectrum (a) and rj distribution (b) for the lepton pair. AlpGen is 
compared with PYTHIA, with and without ME corrections. Relative difference with 
respect to ME corrected PYTHIA is shown for each plot. 



without ME corrections has the shower starting scale set to ^/s. For AlpGen we 
used a minimum for the additional ME generated parton of 20 GeV, and the 
minimum px for the cone algorithm that steers the matching was set to 25 GeV, as 
recommended by the authors. 

AlpGen lepton pair px spectrum appears to be softer than ME-corrected PYTHIA. 
This also translates into a broader t] distribution for AlpGen. If we compare AlpGen 
to both ME corrected and to uncorrected PYTHIA, it appears that AlpGen follows 
uncorrected PYTHIA for low px values and then agrees with corrected PYTHIA in the 
high pt tail. 

We think that this effect is not related to AlpGen itself but rather to PYTHIA. 
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Figure 14: pf spectrum of the electron from Z decay. AlpGen is compared with 
PYTHIA, with and without ME corrections. Relative differences are with respect to 
ME-corrected PYTHIA. 



When PYTHIA is used to shower AlpGen events PYTHIA's native ME corrections are 
switched off, because AlpGen is going to supply such corrections. As noticed in 
Fig. El the side effect of switching off ME corrections in PYTHIA is that the low 
Pt shape changes. AlpGen cannot modify the low pt region, which is completely 
determined by the PS; it can only add high pjn radiation which is actually what it 
does. 

Also the Pt spectrum for the leptons from the Z boson decay (Fig. [HI) shows 
the same behavior. 

The dependency on the resolution cut that separates the ME from the PS re- 
gion is very limited. The lepton pair p^ and rj spectra and the px spectrum for 
the electron from Z are shown in Fig. [TH] for three different choices of the Matrix 
Element cutoff in AlpGen (and correspondingly of the minimum pt of the internal 
cone algorithm): qcut=W, 20, 40 GeV. 

Concerning QCD observables, differential jet rate plots appear to confirm the 
trend observed for letponic observables. Fig. (TH] shows the 1^0 differential jet 
rate. Both ME corrected and uncorrected PYTHIA are shown as a reference. We 
see that AlpGen closely follows uncorrected PYTHIA in the low Qi^o region, then it 
starts to agree with corrected PYTHIA for high values of Qi^o- 

Fig. [T7| shows the 1 — differential jet rate in AlpGen for three different values 
of the ME cutoff in AlpGen. Also ME-corrected PYTHIA is shown as a reference. 
The dependency on the cut is very limited, as already observed for the leptonic 
observables. 
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Figure 15: (a) pt spectrum and (b) r] distribution for the lepton pair, (c) px 
spectrum for the electron. AlpGen has been run with three different values of the 
Matrix Element cutoff. Also ME-corrected PYTHIA is shown as a reference. Relative 
differences are calculated with respect to ME-corrected PYTHIA. 



5.3 AlpGen plus HERWIG 

We tried to shower AlpGen events also with HERWIG. We observed that using HERWIG 
gives results much more consistent with the PYTHIA benchmark. The Z boson trans- 
verse momentum as obtained in AlpGen+HERWIG is shown in Fig. [181 ME corrected 
PYTHIA spectrum is reported in the same plot as a reference. The agreement is 
much better than it was using PYTHIA as a parton shower (Fig. [12] (a)). The low 
Pt region, that was not reproduced in AlpGen+PYTHIA is now well reproduced 
within few percent, except for the very first bins that may be affected by different 
primordial k± tunings. 



6 Conclusion 

We studied the effect of ME corrections in PYTHIA and HERWIG. Both programs can 
fully take into account ME corrections for one additional parton emission. The 
implementation is slightly different. 
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Figure 16: Differential jet rate for the 1^0 transition. AlpGen has been run with 
a ME cutoff of 20 GeV. Both ME corrected and uncorrected (starting shower scale 
set to ^/s) PYTHIA are shown as a reference. Relative differences are with respect 
to ME corrected PYTHIA. 
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Figure 17: Differential jet rates for the transitions 1^0. AlpGen has been run with 
three different values of the Matrix Element cutoff. Also ME-corrected PYTHIA is 
shown as a reference. Relative differences are with respect to ME corrected PYTHIA. 



• PYTHIA modifies the shower in two steps: first the starting scale is raised so 
that any hard emission from the shower is kinematically possible; then the 
emission probability for the first emission is modified to include ME correction 
effect. 

• In HERWIG the shower leaves the phase space for hard emission completely 
uncovered. For this reason the correction is performed in two steps: in the 
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Figure 18: Z boson transverse momentum as obtained in AlpGen+HERWIG and in 
PYTHIA. 

region already filled by the shower the same approach as PYTHIA is used. 
In the remaining part of the phase space the ME for one additional parton 
emission is directly used. 

Both approaches give similar results. When ME corrections are switched off the Z 
Pt spectrum in PYTHIA changes also at low pt, while the low px shape in HERWIG 
remains unchanged. 

We used ME corrected PYTHIA as a test bed for SHERPA and AlpGen. When those 
programs are allowed to emit at most one parton from the ME calculation they 
should give results similar to PYTHIA. Actually SHERPA follows PYTHIA quite well, 
both on lepton and jet observables. On the other hand AlpGen+PYTHIA appears 
to follow uncorrected PYTHIA at low pt, for example it shows softer Z px spectrum 
with respect to ME corrected PYTHIA. This can be traced down to the fact that 
when PYTHIA is used to shower events produced by AlpGen ME corrections are 
switched off. This is done because AlpGen is going to introduce its own corrections. 
The side effect of this is that the low pt shape of the Z px spectrum changes, and 
AlpGen cannot do anything in that region, which is entirely determined by the PS 
alone. AlpGen can only modify the high pt tail of the distribution. 

Using HERWIG to shower AlpGen events turns out to be in much better agreement 
with the PYTHIA benchmark. When using HERWIG to shower AlpGen events native 
HERWIG ME corrections are switched off (like in PYTHIA) but this does not affect the 
low Pt shape, which remains correct. 

The dependency on the cut used to separate the ME and PS regions is limited 
both in AlpGen and in SHERPA. 
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A Generator versions and parameters 



In this work we used PYTHIA version 6.411, SHERPA version 1.1.2, HERWIG++ version 
2.2.1. AlpGen version 2.13 was used together with PYTHIA 6.411 and HERWIG 6.510. 
We used CTEQ6L[21J parton density functions. 

The simulation of multiple interactions has been switched off for all the gener- 
ators; also the QED radiation off final state leptons has been switched off. 

Different settings for ME corrections in PYTHIA presented were obtained using 
parameter MSTP(68). In HERWIG++ we used parameter Evolver :MECorrMode to 
switch on and off ME corrections. 
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